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Abstract 

As a nonbaryonic explanation of massive compact halo objects, a phenomeno- 
logical model is presented which predicts formation of primordial black holes at a 
desired mass scale. The required feature of initial density fluctuation is realized 
making use of the primordially isocurvature fluctuation generated in an inflation- 
ary universe model with multiple scalar fields. 



1 Introduction 



If overdensity of order of unity exists in the hot early universe, a black hole can be 
formed when the perturbed region enters the cosmological horizon. The primor- 
dial black holes (hereafter PBHs) thus produced was a subject of active research 
decades ago (Zel'dovich & Novikov 1967; Hawking 1971) and various observational 
constraints have been obtained against their mass spectrum — ^with no observa- 
tional evidence of their existence at that time (Novikov et al. 1979). 

Recently, however, several independent projects reported observation of massive 
compact halo objects (MACHOs) through gravitational microlensing (Alcock et al. 
1993; Aubourg et al. 1993). It is estimated that their mass is around 0.01 — O.IMq 
and that they occupy ~ 20% of the galactic halo mass which makes up about 
0(10^^) of the critical density (Griest et al. 1995). While the primary candidate 
of MACHOs is substellar baryonic objects such as brown dwarfs, it is difficult to 
reconcile such a large amount of these objects with the observed mass function of 
low mass stars (Richer & Fahlman 1992) and with the infrared observation of dwarf 
component (Boughn & Uson 1995), unless the mass function is extrapolated to the 
lower masses in an extremely peculiar manner. Therefore it is also an interesting 
and potentially important theoretical issue to consider nonbaryonic explanation of 
the origin of MACHOs. 

In the present paper we consider the possibihty that MACHOs consist of PBHs 
produced in the early universe and present a simple model which generates a de- 
sired spectrum of primordial density perturbations in the context of inflationary 
cosmology (Guth 1981; Sato 1981; for a review see, e.g. Olive 1990). In the simplest 
models of inflation with one inflaton scalar fleld, the predicted adiabatic density 
fluctuation has an almost scale-invariant spectrum (Hawking 1982; Starobinsky 
1982; Guth & Pi 1982), unless the inflaton has a peculiar potential. Hence they 
do not predict PBH formation in general. In models with multiple scalar flelds, on 
the other hand, not only adiabatic but also isocurvature fluctuations are generated 
during inflation. The latter can be cosmologically important if energy density of its 
carrier becomes significant in a later epoch (Linde 1985; Kofman & Linde 1987). 
Furthermore it is relatively easier to imprint a nontrivial feature on the spectral 
shape of the isocurvature fiuctuations. Making use of this property here we con- 
struct a model which possesses a peak in the spectrum of total density fluctuation 
at the horizon crossing. Then a significant amount of PBHs can be produced 
around the horizon mass scale when the fiuctuation at the horizon crossing be- 
comes maximal. While our goal is to produce PBHs of mass ~ O.IM© with the 
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abundance of Qbr ~ 10^'^ one can easily see that our model can also be applied 
to PBH formation of different masses and abundances as well by choosing different 
values of model parameters. 

The rest of the paper is organized as follows. In §2 we review basics of PBH 
formation and discuss necessary initial condition of density fluctuations to obtain 
adequate PBHs. Then in §3 possibility of generating the necessary fluctuations 
is considered in the context of inflationary cosmology and a model Lagrangian is 
proposed. §§4-6 are devoted to detailed description of the evolution of the universe 
in this model and in §7 constraints on the model parameters are obtained. §8 is 
the conclusion. 



2 Formation processes of the PBHs. 



PBHs are formed if initial density fluctuations grow sufficiently and a high density 
region collapses within its gravitational radius. First let us review its formation 
process. The background spacetime of the early universe dominated by radiation 
is satisfactorily described by the spatially-flat Priedmann universe, 



ds^ = -dt"^ + a^(t) [dr"^ + (de"^ + sin^ Odif'' 



whose expansion rate is given by 

H\t) 



SttG 



(1) 



(2) 



with p{t) being the background energy density and a dot denotes time derivation. 
Following Carr (1975), let us consider a spherically symmetric high density region 
with its initial radius, R{to), larger than the horizon scale ~ to. The assumption 
of spherical symmetry will be justifled below. Then the perturbed region locally 
constitutes a spatially closed Priedmann universe with a metric. 



ds' = -dt" + Ii\t') 



dr"^ 



+ r\d9^ + sm^9dip'') 



K> 0, 



(3) 



Then the Einstein equation reads 
H'\t') E 



ldR\ _87rG , K 



(4) 



there, where p_|_ is the local energy density. One can choose the coordinate so 
that both the background and the perturbed region have the same expansion rate 
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initially at t = t' = to. Then the initial density contrast, Sq, satisfies 

_ P+Q - Po _ K 

° = — ;; — - Ff2^' 

Po -"o-n-o 

where a subscript implies values at t^. It can be shown that the two time variables 
are related by (Harrison 1970) 

3 dt' dt 

The perturbed region will eventually stop expanding at i? = Rc, which is 
obtained from 

= ^Pod + f^V - ^ = H!(l + *) f - (7) 



3 \RJ Rl ' \RcJ Rl 

as 

Rc^^^^Ro-So^'Ro, (8) 

corresponding to the epoch 

t. ^ 'f. (9) 

do 

The perturbed region must be larger than the Jeans scale, Rj, in order to contract 
further against the pressure gradient, while it should be smaller than the horizon 
scale to avoid formation of a separate universe. We thus require 

Rj^Cste < Rc < te, (10) 

or 

Cs< ^Sl < 1, (11) 

where is the sound velocity equal to 1/ y/S in the radiation dominated era. Since 

1 /2 

RqSq /to is time independent, it suffices to calculate the constraint on 5 at a 
specific epoch, say, when the region enters the Hubble radius, 2R — 2t. We find 
that the amplitude should lie in the range 

^ < 5(R^t) < 1. (12) 

The gravitational radius, Rg, of the perturbed region in the beginning of con- 
traction with Rc ~ Cstc is given by 

Rg = 2GM ~ H^Rl - ^ - clRc < Rc- (13) 
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This is somewhat smaller than but it also implies that a black hole will be 
formed soon after the high density region starts contraction. Thus we expect that 
a black hole with a mass around a horizon mass at t = tc will result and it has 
in fact been shown by numerical calculations (Nadezhin et al. 1978; Bicknell et al. 
1979) that the final mass of the black hole is about (9(10^°'^) times the horizon 
mass at that time. It has been discussed that these black holes do not accrete 
surrounding matter very much and that their mass do not increase even one order 
of magnitude (Carr 1975). Note also that evaporation due to the Hawking radiation 
is unimportant for M > lO^^g (Hawking 1974). 
Since the horizon mass at the time t is given by 

Mhor = 10^ (^) Mq, (14) 



Isec 

what is required in order to produce a significant number of PBHs with mass 
M ~ 0.1 Mq is the sufficient amplitude of density fluctuations on the horizon scale 
at t ~ lO^^sec. Because the initial mass fraction of PBHs, (3, is related with the 
present fraction f^BH as 

/3 = "^^""'T'^ ^bh = 10-«J1bh, (15) 

where teq ~ 10^°sec is the equality time, only an extremely tiny fraction of the 
universe, /3 ~ 10~^^ should collapse into black holes. 

That is, the probability of having a density contrast l/3< S < Ion the 
horizon scale at t ~ 10~^sec should be equal to 0. Let us assume density fluc- 
tuations on the relevant scale obey the Gaussian statistics with the dispersion 
^BH ^ Ij which would be the case in the model introduced in the next sections. 
Then the probability of PBH formation is estimated as 




~ (^BHexp ^- _^ ) , (16) 
which implies that we should have 

5bh ^ 0.05, (17) 
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to produce appropriate amount of PBHs. See figure 1. Although it is true that 
in principle an exponential accuracy is required on the amplitude of fluctuations 
in order to produce the desired amount of PBHs, we have not been able to obtain 
the correspondence between 5bh and Qbk with such an accuracy because numer- 
ical coefficients appearing in the above expressions, such as 18 in (|1^), have been 
calculated based on a rather qualitative argument. We therefore will not attempt 
exceedingly quantitative analysis in what follows. 

Note also that for (5bh = 0.05 the threshold of PBH formation, 6 = 1/3, cor- 
responds to 6.4 standard deviation. It has been argued by Doroshkevich (1970) 
that such a high peak has very likely a spherically symmetric shape. Thus the 
assumption of spherical symmetry in the above discussion is justified and it is also 
expected that gravitational wave produced during PBH formation is negligibly 
small. 



3 Non-flat perturbation in inflationary cosmol- 
ogy 

Since the amplitude of density perturbations on large scales probed by the anisotropy 
of the background radiation (Smoot et al. 1992) is known to be 5 ~ 10^^, the pri- 
mordial fluctuations must have such a spectral shape that it has an amplitude of 
10~^ on large scales, sharply increases by a factor of 10^ on the mass scale of PBHs, 
and decreases again on smaller scales at the time of horizon crossing. It is difficult 
to produce such a spectrum of fluctuations in inflationary cosmology with a single 
component. 

In generic inflationary models with a single scalar field 0, which drives inflation 
with a potential ^ [</>], the root-mean-square amplitude of adiabatic fluctuations 
generated is given by 

If 



on the comoving scale r(</)) when that scale reenters the Hubble radius (Hawking 
1982; Starobinsky 1982; Guth & Pi 1982). The right-hand-side is evaluated when 
the same scale leaves the horizon during inflation. Because of the slow variation 
of and rapid cosmic expansion during inflation, ([T8|) implies an almost scale- 
invariant spectrum in general. Nonetheless one could in principle obtain various 
shapes of fluctuation spectra making use of the nontrivial dependence of 6 (r (</>)) 
on V[(f)] (Hodges & Blumenthal 1990). In order to obtain a desirable spectrum for 
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PBH formation with a mountain on a particular scale, we must employ a scalar 
potential with two breaks and a plateau in between (Ivanov et al. 1994). Such a 
solution is not aesthetically appealing. 

Here we instead consider an inflation model with multiple scalar fields in which 
not only adiabatic but also primordially isocurvature fluctuations are produced. In 
fact it is much easier to imprint nontrivial structure on the isocurvature spectrum 
as mentioned in the beginning. 

We introduce three scalar fields 0i, 02, and 03 in order to generate the desired 
spectrum of density fluctuation. 0i is the inflaton field which induces the new 
inflation (Linde 1982; Albrecht & Steinhardt 1982) with a double- well potential, 
starting its evolution near the origin where its potential is approximated as 

U[^i] ^Vo- ^0t. (19) 

See Linde (1994) and Vilenkin (1994) for the natural realization of its initial con- 
dition. The Hubble parameter during inflation, Hj, is given by 

On the other hand, 02 is a long-lived scalar field which induces primordially isocur- 
vature fluctuations that contribute to black hole formation later. Finally 03 is an 
auxiliary fleld coupled to both 0i and 02 and it changes the effective mass of the 
latter to imprint a speciflc feature on the spectrum of its initial fluctuations. 
We adopt the following model Lagrangian. 

^ = -\{d<t>if - \idhf - \{d<t>zf - V[<t>,, 02, 03] + Ant, (20) 

where £int represents interaction of 0j's with other flelds. Here \^[0i, 02, 03] is the 
effective scalar potential governing the dynamics of the flelds, 

y[01, 02, 03] = t/[0l] + ^(0? - 4>IM + ^03 - ^0302 + ^02 + \^W2, (21) 

where A^, e, 0ic, and 7712 are positive constants. m2 is assumed to be much 
smaller than the scale of inflation, if/, and it does not affect the dynamics of 02 
during inflation. Hence we ignore it for the moment. 

Let us briefly outline how the system evolves before presenting its detailed 
description. In the early inflationary stage 0i is smaller than 0ic, and 03 has 
its potential minimum off the origin. Then 02 also settles down to a nontrivial 
minimum, where it can have an effective mass larger than Hj so that its quantum 
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fluctuation is suppressed. As 0i becomes larger than 0ic, 03 rolls down to its 
origin. Then the potential of 02 also becomes convex and its amplitude gradually 
decreases due to its quartic term. However, since its potential is now nearly flat, 
its motion is extremely slow with its effective mass smaller than Hj well until 
the end of inflation. In this stage quantum fluctuations are generated to 02 with 
a nearly scale-invariant spectrum. Thus the initial spectrum of the isocurvature 
fluctuations due to 02 has a scale- invariant spectrum with a cut-off on a large scale. 

After inflation, the Hubble parameter starts to decrease in the reheating pro- 
cesses. As it becomes smaller than the effective mass of 02, the latter starts rapid 
coherent oscillation. 02 dissipates its energy in the same way as radiation in the 
beginning when its oscillation is governed by the quartic term. But later on when 
A20| becomes smaller than m|, its energy density decreases more slowly in the 
same manner as nonrelativistic matter. Thus 02 contributes to the total energy 
density more and more later, which implies that the total density fluctuation due 
to the primordially isocurvature fluctuations or 02 grows with time. Since what is 
relevant for PBH formation is the magnitude of fluctuations at the horizon cross- 
ing, we thus obtain a spectrum with a larger amplitude on larger scale until the 
cut-off scale in the initial spectrum is reached, that is, it has a single peak on the 
mass scale of PBH formation. 

In the above scenario we have assumed that 02 survives until after the PBH 
formation. On the other hand, were 02 stable, it would soon dominate the total 
energy density of the universe in conflict with the successful nucleosynthesis. As 
a natural possibility we assume that 02 decays through gravitational interaction, 
so that it does not leave any unwanted relics with its only trace being the tiny 
amount of PBHs produced. 

In the subsequent sections we describe the detailed evolution of the above model 
and obtain constraints on the model parameters to produce the right amount of 
PBHs on the right scale. 

4 Background evolution 

First we consider the evolution of the homogeneous part of the flelds. During 
inflation, the behavior of the inflaton is governed by the C/[0i] part of the potential. 
Solving the equation of motion with the slow-roll approximation. 




(22) 
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we find 

= T' (23) 

where (pu is the field amphtude at some initial epoch ti. The above approximate 
solution remains valid until |f/"[0i]| becomes as large as 9if| at t = tj, when 
infiationary expansion is terminated and we find (pi{tf) = 3Hj/Xi. Then (^) can 
also be written as 

i2f \ 3Hj 3H^ 3Hj 

where r(t) is the e-folding number of exponential expansion after t (< t/) and the 
last approximation is valid when r ^ 1. From now on, we often use r(t) as a new 
time variable or to refer to the comoving scale leaving the Hubble radius at t. Note 
that it is a decreasing function of t. 

As stated in the last section, we are taking a view that (pi determines fate of 
03 and that 03 controls evolution of 02 but not vice versa. In order that 03 does 
not affect evolution of 0i the inequality 

A1A3 > e\ (25) 

must be satisfied, while we must have 

A2A3 > (26) 

so that 02 does not affect the motion of 03. We assume these inequalities hold 
below. 

When 01 < 01c, both 02 and 03 have nontrivial minima which we denote by 
02m and 03m, respectively. From 

V2 = -1^0302 + A20^ = 0, (27) 

and 

Vs = e(0? - 02J03 + A30^ - i/0^03 = 0, (28) 
with Vj = dV/d(j)j, we find 

A20L(t) = ^0L(i) (29) 
A30LW = e(0L-0?(^))+^02m(t) = e(0L-0?(t)), (30) 



9 



where ( PB] ) was used in the last expression. In the early inflationary stage when 
(f>i -C 01c, the effective mass-squared of (pj, Vjj, at the potential minimum is given 
by 



2A202m = ^Wlc 
^3 



33 KA'l, <A'2m, <^3m 



3e 



2z/e 



3z/e 



where Tc is the epoch when 



We choose parameters such that 



(31) 
(32) 



ue > AiAsTc, and e > Ait^ 



(33) 



Then V22 and V33 are larger than Hj initially at the potential minimum, so that 
both 02 and ^3 settle down to 02m(^) and (t>3mit), respectively. 

03m(^) decreases down to zero at r = Tc when V33[03m] also vanishes. Then 
^3 [03m = 0] starts to increase according to 0i and soon acquires a large positive 
value, which implies that 03 practically traces the evolution of 03m(i) down to 
zero without delay. On the other hand, 02 evolves somewhat differently because 
it does not acquire a positive effective mass from 03 at the origin. Although 02 (t) 
traces 02m (i^^) initially, as V22[02m] becomes smaller it can no longer catch up with 
02m (^)- From a generic property of a scaler field with a small mass in the De Sitter 
background, one can show that this happens when the inequality 



d(f)2m{t) 



02m (t) dt 



> 



gets satisfied, or at 



with 



r 



(Plin) = 0i 



'A1A3 

2z/e 



22 [V2m\ 

3H 



n = Tc 



ve 



(34) 



(35) 



(36) 



y 2A1A3 2A2r; 

Thus 02 slows down its evolution. In the meantime 03 vanishes. Then 02 is 
governed by the quartic potential. We can therefore summarize its evolution during 
inflation as 



2A1A2A3 



3Hj 



in-r) 



1 -1 



< r < n. 



(37) 
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Since A202/ is adequately smaller than 02 remains practically constant in the 
latter regime. Let us also write down time dependence of its effective mass-squared 
for later use. 

ueHj ( 1 



A1A3 Vt-c ' ' ~ 



'22\ 



3A20I 



l + ^(r.-r)l" 



(3^ 



< r < Ti. 



5 Generation of fluctuations 



We now consider fluctuations in both the scalar fields and metric variables in a 
consistent manner. We adopt Bardeen's (1980) gauge-invariant variables (^a and 
with which the perturbed metric can be written as 

ds^ = -(1 + 2<!?A)dt^ + a{tf{l + 2^H)dx^, (39) 

in the longitudinal gauge. In this gauge scalar field fluctuation 6(j)j coincides with 
the corresponding gauge-invariant variable by itself. 

Assuming an exp(ikx) spatial dependence and working in the Fourier space, 
the perturbed Einstein and scaler field equations are given by 

^A + ^H = (40) 

^H + H^H = -47rG'(0i(50i + 02(502 + 03^03) (41) 



b, + 3H64>, + ( ^ + ^i^) ^"^J = + *A0, - 3<1''^^0, - J2 "^J^^^i^ (42) 



Note that all the fluctuation variables are functions of k and t. 

The above system is quite complicated at a glance. However, using constraints 
on various model parameters we have obtained so far, i.e. (pSf), (pBf), and (p3D, 
it can somewhat be simplified. First, since 03 has an effective mass larger than 
Hj during inflation except in the vicinity of r = Tc, quantum fluctuation on 503 
is suppressed and moreover its energy density practically vanishes by the end of 
inflation. We can therefore neglect fluctuations in 03. On the other hand, we can 
show that 

' ' (43) 



V A2A3 

with the help of (p6D and (^). Hence (|4l|) and (^2|) with j = 1 reduce to 

+ H^H = -47rG0i50i (44) 
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5^1 + 3HS'(Pi + l^-^ + Vuj Hi = 2Vi^A - 4$/^0i. (45) 

Thus only (pi contributes to adiabatic fluctuations and it can be calculated in the 
same manner as in the new inflation model with a single scalar fleld. This is as 
expected because 0i dominates the energy density during inflation. In fact since 
we are only interested in the growing mode on the super-horizon regime which 
turns out to be weakly time-dependent as can be seen from the flnal result, we can 
consistently neglect time derivatives of metric perturbations and terms with two 
time derivatives in and (^) during inflation. We thus flnd 



= -^A = (46) 

The resultant amplitude of scale-invariant adiabatic fluctuations depends on Ai 
and one can normalize its value using the COBE observation (Smoot et al. 1992) 
as 

Al = 1.3 X 10"^^ (47) 

On the other hand, 502 satisfies (^) with j = 2. From quantum field theory in 
De Sitter spacetime, it has a root-mean-square amplitude 502 = {H^/2k^Yl^ when 
the /c-mode leaves the Hubble radius if V22 is not too large. Since V2 vanishes when 
02 = 02m and V2$A remains small even for r < r/, we can neglect all the terms in 
the right-hand side, to yield 

502 + 3if7502 + ^22502 = 0, (48) 

when k <^ a{t)Hj. We can find a WKB solution with the appropriate initial 
condition, 




where tk is the time when fc-mode leaves the Hubble radius: k = a{tk)Hi. The 
above expression is valid when \S\ <C S"^. In terms of r, (ESI) can be expressed as 



IMt, r)^^§ ' exp {/; S(ryir' -'-(r.-r)}, (50) 
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S{r) 



t 1 



1 

-1\ 2 



< r < Tu 



where = ritk) and r^, = ^^^^^ . The above equahty is vahd until the end of 
inflation at t = t/ or r = 0. 



6 Evolution of the universe after inflation 

Let us assume the universe is rapidly and efficiently reheated att = tf for simplicity 
to avoid further complexity, (see, e.g. Kofman et al. (1994), Shtanov et al. (1995), 
and Boyanovsky et al. (1995) for recent discussion on efficient reheating.) Then 
the reheat temperature is given by 



Tr = O.l^HjMpi. (51) 

If there is no further significant entropy production later, one can calculate the 
epoch, r(L), when the comoving length scale corresponding to L pc today left the 
Hubble horizon during inflation as 

r(L) = 37 + In ( + 1 In ( . (52) 



Ipc; 2 VlQiOGeV 

Then the comoving horizon scale at t = lO^^sec, or L = 0.03 pc corresponds to 
r = 34 = Tin and the present horizon scale ^ 3000Mpc to r = 59. 
On the other hand, 02 (^) and 502 (k,t) evolve according to 

02 + 3if 02 + A202 + m^02 = 0, (53) 
(502 + 3i/502 + (3A20^ + m2)502 ^ 0, (54) 

where the Hubble parameter is now time-dependent: H = l/2t, and the latter 
equation is valid for k <^ aH. 

When becomes smaller than A202 (^ ^2)5 both 02 and 502 start rapid 
oscillation around the origin. Using (|50D one can express the amplitude of the 
gauge- invariant comoving fractional density perturbation of 02, A2, as 



A2 = — ^02502 - 02(502 - 02^"^"^) ^ 4^ 
P2 V / 02 



T=0 



.4^, (55) 

02f 
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in the beginning of oscillation. Here 

P2 = '^<P2 + + ^"^202 (56) 

is the energy density of 02- 

Using the virial theorem one can easily show that it decreases in proportion to 
a~^{t) as long as A202 ^ ^'i- Thus the amplitude of 02 decreases with a~^{t). 
On the other hand, 502 has a rapidly oscillating mass term when A202 > 
which causes parametric amplification. We have numerically solved equations ( |53D 
and (Q) with various initial conditions with |02| ^ |'^02| initially. We have found 
that in all cases the amplitude of 502 remains constant as long as m2 is negligible. 
Thus A2 increases in proportion to t^^^ in this regime and A2 becomes as large as 



k2 



A,^4^^ ^) , ,57) 

02f \ mi ) 

while the ratio of p2 to the total energy density, ptot, which is now dominated by 
radiation, remains constant: 

Ptot \Mpi) 

As A202 becomes smaller than 02 and 502 come to satisfy the same equation 
of motion, see (|53D and (^41), and A2 saturates to the constant value (|57D. At the 
same time p2 starts to decrease less rapidly than radiation, in proportion to a^^{t). 
Since A2 contributes to the total comoving density fluctuation by the amplitude 

A = -^A2, (59) 

Ptot 

it increases in proportion to a{t) oc t^^^. In the beginning of this stage, we find 
^ m^/A20^f, to yield 




A^8n-^l-^\ i^P-] , (60) 



at a later time t. 

In order to relate it with the initial condition required for PBH formation, we 
must estimate it at the time fc-mode reenters the Hubble radius, tl, defined by 

k = 2na{tl)H{tl) = (61) 
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Since k can also be expressed as k = 2'7ra/e ^^Hj, the amplitude of comoving 
density fluctuation at t = is given by 

7 Constraints on model parameters 

In §2 we discussed the necessary condition on the amplitude of fluctuations for 
PBH formation using the uniform Hubble constant gauge. Hence we should calcu- 
late the predicted amplitude in this gauge, which is a linear combination of A and 
the gauge-invariant velocity perturbation. However, in the present case in which 
A grows in proportion to a{t) in the radiation-dominant universe, one finds that 
the latter quantity vanishes and that density fluctuation in the uniform Hubble 
constant gauge coincides with A (Kodama & Sasaki 1984). Thus we finally ob- 
tain the quantity to be compared with 6bh in ([T^), namely, the root-mean-square 
amplitude of density fluctuation on scale r = 27i/k at the horizon crossing, S{r), 
as 



6{r)- 

with 



2f 



(27r) 

[ sio) ) 

where we have used (150|). We also find 



e"^C/(rk,re,r,), (63) 



^1 





(65) 



from (|3B|) and (0). 

The remaining task is to choose values of parameters so that 6{r) has a peak on 
the comoving horizon scale at t = 10~^sec, which we denote by rm, corresponding 
to Tk = Tm = 34, with its amplitude 5{rm) = 0.05. We thus require 

dln6{r) _ S'{r^) 1 

+ + -- (66) 



4rk [TcTk (tc - Tk)] 2V Tc Tk/ 2 
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vanishes at = Tm, which gives us a relation between and r^,. Since roughly 
corresponds to the comoving scale where scale-invariance of primordial fluctuation 
A2 is broken, the peak at 5(rm) becomes the sharper, the closer Tc approaches Tm- 
For example, if we take = 30 we find 

4i/e 



3Ai A2 



200, (67) 



so that 



Jj2 

C/ = 0.13 and = 0.045— (68) 

A2 



In order to have 6{7\^i) = 0.05, we find 

1 ^ _ 1 ^ in-l5 



\Mpi 



= 1.7 X 10-'^ (69) 



which can easily be satisfied with some reasonable choices of A2 and Hj. However, 
it is not the final constraint. Since we are assuming that the universe is dominated 
by radiation at this time, we require 

^ = 2vr f y f ^i^") ' e- < 1. (70) 

Furthermore 02 should decay some time after t — 10~^sec so as not to dominate the 
energy density of the universe which would hamper the primordial nucleosynthesis. 
Assuming that it decays only through gravitational interaction, its life time is given 

-*.-# = l<'-"fl5^)"-. ,71) 



m 



2 



Now we have displayed all the necessary equalities and inequalities the model 
parameters should satisfy. Since there is a wide range of allowed region in the multi- 
dimensional space of parameters, we do not work out the details of the constraints 
but simply give one example of their values with which all the requirements are 
satisfied: 



Hi 


= 1.7 X lO^^GeV, 


1712 


= 3.2 X lO^GeV, 


Al 


= 1.3 X 10-^^ 


A2 


= 1.4 X 10-^ 


A3 


= = 6.7 X 10-^ 


e 


= 2.0 X 10-^\ 



(72) 
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for which P2/ Ptot = 0.1 at t = 10~®sec and inequahties (P3| ) and (^) are maximally 
satisfied. 

In figure 2 we have depicted the qualitative mass spectrum of produced PBHs 
for different values of where 

/3(M) = 5(r)exp(-^^) (73) 

has been shown as a function of the horizon mass when the scale r reenters the 
horizon. 

8 Conclusion 

In the present paper we have considered possibility to produce a significant amount 
of PBHs on a specific mass scale by generating appropriate spectrum of density 
fiuctuations in inflationary cosmology. We have reached a model with the desired 
feature making use of a simple polynomial potential (pTD without introducing any 
break in the potential of the scalar fields. We have chosen values of the model 
parameters so that these PBHs can account for the observed MACHOs. In order 
to set the order of magnitude of the mass scale of the black holes and that of their 
abundance correctly, we had to tune some combinations of model parameters such 
as (|67|) and (|69|) with two digits' accuracy. However, there exists a wide range of 
allowed region in the parameter space to realize it. We also note that the precise 
values such as those quoted in (^) are not of much significance, primarily because 
the formula for the fraction of PBHs (0) is only a qualitative one. In this respect 
we have restricted ourselves to the analytic treatment of evolution of fiuctuation, 
which is not exponentially accurate. In the event a more precise formula for PBH 
fraction is obtained, full numerical analysis of fiuctuations would be required. At 
the present stage, however, analytic treatment is more appropriate with which 
dependence of the results on physical parameters are understood more clearly. 

It is evident that our model can be applicable to produce PBHs with a different 
mass and abundance by slightly changing values of parameters. For example, we 
could produce black holes with mass ~ lO^M© which would act as a central engine 
of AGNs. These black holes are usually considered to have formed in the post- 
recombination universe (Loeb 1993), but they might have formed in the early 
universe at t ~ lOsec corresponding to the onset of primordial nucleosynthesis. 
Note that this would not hamper successful nucleosynthesis because the root-mean- 
square amplitude of density fluctuation required for such black hole formation is 



17 



still much smaller than unity. 
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Figure captions 

Figure 1 Fraction of primordial black holes as a function of root-mean-square am- 
plitude of density fluctuations (eq.JlBl) Gaussian distribution of fluctuations 
is assumed. 

Figure 2 Expected mass spectrum of primordial black holes with different values 

of Tc. 
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